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The complete genomes of Thermus oshimai JL-2 and T. thermophilus JL-18 each consist of a cir- 
cular chromosome, 2.07 Mb and 1.9 Mb, respectively, and two plasmids ranging from 0.27 Mb 
to 57.2 kb. Comparison of the T. thermophilus JL-18 chromosome with those from other strains of 
T. thermophilus revealed a high degree of synteny, whereas the megaplasmids from the same 
strains were highly plastic. The T. oshimai JL-2 chromosome and megaplasmids shared little or no 
synteny with other sequenced Thermus strains. Phylogenomic analyses using a concatenated set 
of conserved proteins confirmed the phylogenetic and taxonomic assignments based on 1 6S 
rRNA phylogenetics. Both chromosomes encode a complete glycolysis, tricarboxylic acid (TCA) 
cycle, and pentose phosphate pathway plus glucosidases, glycosidases, proteases, and peptidases, 
highlighting highly versatile heterotrophic capabilities. Megaplasmids of both strains contained a 
gene cluster encoding enzymes predicted to catalyze the sequential reduction of nitrate to nitrous 
oxide; however, the nitrous oxide reductase required for the terminal step in denitrification was 
absent, consistent with their incomplete denitrification phenotypes. A sox gene cluster was identi- 
fied in both chromosomes, suggesting a mode of chemolithotrophy. In addition, n/f and psr gene 
clusters in T. oshmai JL-2 suggest respiratory nitrite ammonification and polysulfide reduction as 
possible modes of anaerobic respiration. 

Abbreviations: NCBI- National Center for Biotechnology Information (Bethesda, MD, USA), 
IMG- JGI Integrated Microbial Resource 



Introduction 



The Great Boiling Spring (GBS] geothermal system 
is located in the northwestern Great Basin near the 
town of Gerlach, Nevada. Geothermal activity is 
driven by deep circulation of meteoric water, which 
rises along range-front faults at temperatures up to 
96 Q C. A considerable volume of geomicrobiology 
research has been conducted in the GBS system, 
including coordinated cultivation-independent mi- 
crobiology and geochemistry studies [1-4], habitat 
niche modeling [3], thermodynamic modeling [1,5], 
microbial cultivation and physiology [6,7], and in- 
tegrated studies of the nitrogen biogeochemical 
cycle (N-cycle [5,6,8]]. The latter group of studies is 
arguably the most detailed body of work on the N- 
cycle in any geothermal system. Those studies re- 
vealed a dissimilatory N-cycle based on oxidation 



In high temperature sources such as GBS and 
Sandy's Spring West (SSW], ammonia oxidation 
occurs at temperatures up to at least 82 Q C at rates 
comparable to those in nonthermal aquatic sedi- 
ments [5]. Several lines of evidence, including deep 
16S rRNA gene pyrosequencing datasets and quan- 
titative PCR, suggest ammonia oxidation is carried 
out by a single species of ammonia-oxidizing 
archaea closely related to "Candidatus 
Nitrosocaldus yellowstonii", which comprises a 
substantial proportion of the sediment microbial 
community in some parts of the springs [5,9]. Ni- 
trite oxidation appears to be sluggish or non- 
existent in the high temperature source pools since 



and subsequent denitrification of ammonia sup- 
plied in the geothermal source water. 
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nitrite accumulates in these systems and 16S rRNA 
gene sequences for nitrite-oxidizing bacteria have 
not been detected in clone library and pyrotag cen- 
suses [1,5]. Finally, the nitrite and nitrate that are 
produced are denitrified in the sediments to both 
nitrous oxide and dinitrogen; however, a high flux 
of nitrous oxide, particularly in the ~80 Q C source 
pool of GBS, suggested the importance of incom- 
plete denitrifiers [6] and electron donor stimula- 
tion experiments suggested a key role for hetero- 
trophic denitrifiers [5]. 

A subsequent cultivation study of heterotrophic 
denitrifiers in GBS and SSW resulted in the isola- 
tion of a large number of denitrifiers belonging to 
Thermus thermophilus and T. oshimai, including 
strains T. oshimai JL-2 and T. thermophilus JL-18 
[6]. Strikingly, although Thermus strains were iso- 
lated using four different isolation strategies, nine 
different electron donor/acceptor combinations, 
and four different sampling dates, all isolates of 
these two species were able to convert nitrate-N 
stoichiometrically to nitrous oxide-N, but appeared 
unable to reduce nitrous oxide to dinitrogen. This 
physiology, combined with high nitrous oxide flux- 
es in situ suggested a significant role of T. oshimai 
and T. thermophilus in the unusual N-cycle in these 



hot springs. However, the genetic basis of this phe- 
notype remained unknown. Here we present the 
complete genome sequences of T. oshimai JL-2 and 
T. thermophilus JL-18, compare them to genomes of 
other sequenced Thermus spp., and discuss them 
within the context of their potential impacts on bi- 
ogeochemical cycling of carbon, nitrogen, sulfur, 
and iron. 

Classification and features 

The genus Thermus currently comprises 16 species 
and includes the well-known T. aquaticus and the 
genetically tractable T. thermophilus. The genome 
of T. oshimai JL-2 is the first finished genome to be 
reported from that species, while T. thermophilus 
JL-18 is the fourth genome to be sequenced from 
that species, the other being T. thermophilus HB27, 
HB8, and SG0.5JP17-16. Figure 1 shows the rela- 
tionship of T. oshimai JL-2 and T. thermophilus JL- 
18 to other Thermus species, as determined by 
phylogenomic analysis of highly conserved genes, 
which supports the taxonomic identities previously 
determined by 16S rRNA gene phylogenetic analy- 
sis [6]. Table 1 shows general features of T. oshimai 
JL-2 and T. thermophilus JL-18. 
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Figure 1. Phylogenomic tree highlighting the position of Thermus 
oshimai JL-2 and Thermus thermophilus JL-18. Thirty-one bacterial phy- 
logenetic markers were identified using Amphora [10]. Maximum- 
likelihood analysis was carried out with a concatenated alignment of all 
31 proteins using RAxML Version 7.2.6 [1 1] and the tree was visualized 
using iTOL [12]. Red circles indicate bootstrap support >80% (100 rep- 
licates). Scale bar indicates 0.1 substitutions per position. The protein 
FASTA files for all the species are from NCBI, except for the following 
species, which are from IMG: Thermus igniterrae ATCC 700962 (Taxon 
OID: 2515935625), Thermus oshimai DSM 12092 (Taxon OID: 
2515463139), Thermus oshimai JL-2 (Taxon OID: 2508706991), 
Thermus sp. RLM (Taxon OID: 2514335427). 
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Table 1(a). Classification and general features of Thermus oshimai JL-2 according to the MIGS recommendations [13]. 



MIGS ID 


Property 


Term 




Evidence code 3 










TAS [14] 






Domain Bacteria 












TAS [15] 






Phylum Deinococcus-Thermus 








Class Deinococci 




TAS [16,17] 










. .,. 
Current classification 






TAS [16,18] 






Order Thermales 












TAS [16,19] 






Family Thermaceae 












TAS [20-22] 






Genus Thermus 












TAS [23] 






Species Thermus oshimai 












TAS [6] 






Type strain JL-2 






Gram stain 






tac; n ~>.) 






Negative 






Cell shape 


Rod 




1 r\o [0,z.s> J 




Motility 


Non-moti le 




NAS [13] 








Sporulation 






TAS [13] 




Nonsporulating 








Temperature range 


Not reported 














Optimum temperature 


70 °C 




TAS [13] 




Carbon source 


oeverai mono- anu uisaccnariues, 


some organic 


TAS [13] 










acids and amino acids 






Energy source 






TAS [6,23] 






Chemoorganotroph 






Terminal electron acceptor 


0 2 , NO," 




TAS [6,23] 










naDiiai 


Terrestrial hot springs 




TAS [6,23] 


MIGS-6 








Salinity 






Tic n l 
i ad [ i j 


MIGS-6.3 


3.90 g/L total dissolved solids 






Oxygen 






TAS [6,23] 


MIGS-22 


Facultative anaerobe (nitrate reduction) 




Biotic relationship 






TAS [6,23] 


MIGS-15 




Free living 






Pathogenicity 






NAS 


MIGS-14 


Non-pathogenic 






MIGS-4 


Geographic location 


Sandy's Spring West, Great Boi 1 inj 


I Springs geo- 


TAS [6] 




thermal field, Nevada 








Sample collection time 


October, 2008 




TAS [6] 


MIGS-5 






Latitude 






1 Ao [ 1 J 


MIGS-4. 1 




N40° 39.182' 






Longitude 








MIGS-4.2 




W119° 22.496' 








Depth 


Sediment/water interface (shallow^ 




TAS [1] 


MIGS-4. 3 


) 




Altitude 






NAS 


MIGS-4.4 




1,203 m 







""Evidence codes - IDA: Inferred from Direct Assay; TAS: Traceable Author Statement (i.e., a direct report exists in the 
literature); NAS: Non-traceable Author Statement (i.e., not directly observed for the living, isolated sample, but based 
on a generally accepted property for the species, or anecdotal evidence). These evidence codes are from Gene Ontol- 
ogy project [24] . 
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MIGS ID 


Property 


Term 


Evidence code 3 






Domain Bacteria 


TAS [14] 






Phylum Deinococcus-Thermus 


TAS [15] 






Class Deinococci 


TAS [16,17] 




Current classification 


Order Thermales 

Family Thermaceae 

Genus Thermus 

Species Thermus thermophilus 

Type strain JL-1 8 


TAS [16,18] 
TAS [16,19] 
TAS [20-22] 
TAS [25-27] 
TAS [28] 




Gram stain 


Negative 


TAS [28] 




Cell shape 


Rod 


TAS [6,28] 




Motility 


Non-motile 


TAS [28] 




Sporulation 


Nonsporulating 


TAS [28] 




Temperature range 


Not reported 






Onfimiim fpmnerafurp 

V . V IJ 111 1 IUI 1 1 LV, 1 1 1 ij 1 ULUI L . 


70 °C 


TAS [28] 




^aiuuii suurce 


Several mono- and disaccharides; some organic acids 
and amino acids 


TAS [28] 




Energy source 


Chemoorganotroph 


TAS [28] 




Terminal electron acceptor 


0 2 , NO," 


TAS [6] 
TAS [6] 


MICS-6 


Habitat 


Terrestrial hot springs 


MIGS-6.3 


Salinity 


3.90 g/L total dissolved solids 


TAS [1] 


MIGS-22 


Oxygen 


Facultative anaerobe (nitrate reduction) 


TAS [6,13] 


MIGS-15 


Biotic relationship 


Free living 


TAS [6,13] 


MIGS-14 


Pathogenicity 


Non-pathogenic 


NAS 


MIGS-4 


Cipopranhir location 


Sandy's Spring West, Great Boiling Springs 
geothermal field, Nevada 


TAS [6] 


MIGS-5 


Sample collection time 


12/2008 


TAS [6] 


MIGS-4.1 


Latitude 


N40° 39.182' 


TAS [1] 


MIGS-4.2 


Longitude 


W119° 22.506' 


TAS [1] 
NAS 


MIGS-4.3 


Depth 


Sediment/water interface (shallow) 


MIGS-4.4 


Altitude 


1,203 m 



Evidence codes - IDA: Inferred from Direct Assay; TAS: Traceable Author Statement (i.e., a direct report exists in the 
literature); NAS: Non-traceable Author Statement (i.e., not directly observed for the living, isolated sample, but based 
on a generally accepted property for the species, or anecdotal evidence). These evidence codes are from Gene On- 
tology project [24]. 
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Genome sequencing information 

Genome project history 

T. oshimai JL-2 and T. thermophilus JL-18 were se- 
lected based on their important roles in 
denitrification and also for their biotechnological 
potential. The genome projects for both the organ- 
isms are deposited in the Genomes Online Data- 
base [29] and the complete sequences are depos- 
ited in GenBank. Sequencing, finishing, and anno- 
tation were performed by the DOE Joint Genome 
Institute QGI). A summary of the project and in- 
formation associated with MIGS version 2.0 



compliance [13] are shown [T. oshimai JL-2; Table 
2 (a] and T. thermophilus JL-18; Table 2(b]]. 

Growth conditions and DNA isolation 

Axenic cultures of T. oshimai JL-2 and T. 
thermophilus JL-18 were grown aerobically on 
Thermus medium as described [6] and DNA was 
isolated from 0.5-1.0 g of cells using the Joint Ge- 
nome Institute's QGI] cetyltrimethyl ammonium 
bromide protocol [30]. 



Table 2(a). Thermus oshimai JL-2 genome sequencing project information 
MIGS ID Property Term 



MIGS-31 Finishing quality 

MIGS-28 Libraries used 

MIGS-29 Sequencing platforms 

MIGS-31 .2 Fold coverage 

MIGS-30 Assemblers 

MIGS-32 Gene calling method 

Genome Date of Release 



Finished 

454 standard and PE, lllumina 
lllumina GAM, 454-GS-FLX-Titanium 
38.3x (454), 2,228.9x (lllumina) 
Newbler v 2.3 (pre-release) 
Prodigal 1 .4, GenePRIMP 



Genbank ID 



CP003249.1 (chromosome) 
CP003250.1 (Plasmid pTHEOSOl) 
CP003251.1 (Plasmid pTHEOS02) 

Genbank Date of Release November 5, 2012 

GOLD ID Gc02356 

Project relevance Biotechnological 



Table 2(b). Thermus thermophilus JL-18 genome sequencing project information 
MIGS ID Property Term 



MIGS-31 Finishing quality 

MIGS-28 Libraries used 

MIGS-29 Sequencing platforms 

MIGS-31. 2 Fold coverage 

MIGS-30 Assemblers 

MIGS-32 Gene calling method 



Finished 

454 standard and PE, lllumina 
lllumina GAM, 454-GS-FLX-Titanium 
38.1 x (454), 300x (lllumina) 
Newbler v 2.3 (pre-release) 
Prodigal 1.4, GenePRIMP 



Genome Date of Release Oct 21, 2011 



Genbank ID 



CP003252.1 (chromosome) 
CP003253.1 (plasmid pTTJLI 801) 
CP003254.1 (plasmid pTTJLI 802) 

Genbank Date of Release April 9,2012 

GOLD ID Gc02194 

Project relevance Biotechnological 



http://standardsingenomics.org 



453 



Thermus oshimai JL-2 and T. thermophilus JL-18 



Genome sequencing and assembly 

The draft genomes of Thermus oshimai JL-2 and 
Thermus thermophilus JL-18 were generated at the 
DOE Joint Genome Institute QGI] using a combina- 
tion of Illumina [31] and 454 technologies [32]. 

For T. oshimai JL-2, we constructed and sequenced 
an Illumina GAii shotgun library which generated 
146,341,736 reads totaling 11,122 Mb, a 454 Titani- 
um standard library which generated 181,476 reads 
and 1 paired end 454 library with an average insert 
size of 8 kb that generated 285,154 reads totaling 
146.6 Mb of 454 data. For T. thermophilus JL-18, we 
constructed and sequenced an Illumina GAii shotgun 
library that generated 74,093,820 reads totaling 
5,631.1 Mb, a 454 Titanium standard library that 
generated 212,217 reads and 1 paired end 454 li- 
brary with an average insert size of 7 kb that gener- 
ated 121,082 reads totaling 116.9 Mb of 454 data. 
All general aspects of library construction and se- 
quencing performed at the JGI can be found at [30]. 
The initial draft assemblies of T. oshimai JL-2 and T 
thermophilus JL-18 contained 39 contigs in 2 scaf- 
folds and 75 contigs in 3 scaffolds, respectively. 

The 454 Titanium standard data and the 454 paired 
end data were assembled together with Newbler, 
version 2.3-PreRelease-6/30/2009. The Newbler 
consensus sequences were computationally shred- 
ded into 2 kb overlapping fake reads (shreds]. 
Illumina sequencing data was assembled with VEL- 
VET, version 1.0.13 [33], and the consensus se- 
quence were computationally shredded into 1.5 kb 
overlapping fake reads (shreds]. We integrated the 
454 Newbler consensus shreds, the Illumina VEL- 
VET consensus shreds and the read pairs in the 454 
paired end library using parallel phrap, version SPS - 
4.24 (High Performance Software, LLC]. The soft- 
ware Consed [34] was used in the following finishing 
process. Illumina data was used to correct potential 
base errors and increase consensus quality using the 
software Polisher developed at JGI (Alia Lapidus, 
unpublished]. Possible mis-assemblies were cor- 
rected using gapResolution (Cliff Han, unpublished], 
Dupfmisher [35] or sequencing cloned bridging PCR 
fragments with subcloning. Gaps between contigs 
were closed by editing in Consed, by PCR and by 
Bubble PCR (J-F Cheng, unpublished] primer walks. 
Additional reactions were necessary to close gaps 
and to raise the quality of the finished sequence (7! 
oshimai JL-2: 20 reactions; T. thermophilus JL-18: 
45]. 

The total size of the genomes are 2,401,329 bp (7. 
oshimai JL-2) and 2,311,212 bp [T. thermophilus 



JL-18). The final assembly of T. oshimai JL-2 genome 
is based on 91.8 Mb of 454 draft data which pro- 
vides an average 38.3x coverage of the genome and 
5,349.4 Mb of Illumina draft data which provides an 
average 2,228.9x coverage of the genome. The final 
assembly of T. thermophilus JL-18 genome is based 
on 87.7 Mb of 454 draft data which provides an av- 
erage 38.1x coverage of the genome and 690 Mb of 
Illumina draft data which provides an average 300x 
coverage of the genome. The data and metadata are 
made available at the JGI Integrated Microbial Re- 
source website (IMG) [31]. 

Genome annotation 

Initial identification of genes was done using Prodi- 
gal [36], a part of the DOE-JGI Annotation pipeline, 
followed by manual curation using GenePRIMP 
[37]. The predicted ORFs were translated into puta- 
tive protein sequences and searched against data- 
bases including: NCBI nr, Uniprot, TIGR-Fam, Pfam, 
PRIAM, KEGG, COG, and Interpro. Additional anno- 
tations and curations were performed using the 
Integrated Microbial Genomes - Expert Review 
(IMG-ER) platform [33]. 

Genome properties 

The T. oshimai JL-2 genome includes one circular 
chromosome of 2,072,393 bp (2205 predicted 
genes), a circular megaplasmid, pTHEOSOl (0.27 
Mb, 268 predicted genes), and a smaller circular 
plasmid, pTHEOS02 (57.2 Kb, 75 predicted genes), 
for a total size of 2,401,329 bp. Of the total 2,548 
predicted genes, 2,488 were protein-coding genes. 
A total of 2,015 (79%) protein-coding genes were 
assigned to a putative function with the remaining 
annotated as hypothetical proteins. The properties 
and the statistics of the genome are summarized in 
Table 3a, Table 3b, Table 3c and Figure 2). 

The T. thermophilus JL-18 genome includes one cir- 
cular chromosome of 1,902,595 bp (2,057 predict- 
ed genes), a circular megaplsmid, pTTJL1801 (0.26 
Mb, 279 predicted genes), and a smaller circular 
plasmid, pTTJL1802 (0.14 Mb, 172 predicted 
genes), for a total size of 2,311,212 bp. Of the total 
2,508 predicted genes, 2,452 were protein-coding 
genes. A total of 1,979 (79%) of protein-coding 
genes were assigned to a putative function with the 
remaining annotated as hypothetical proteins. The 
properties and the statistics of the genome are 
summarized in Table 4a, Table 4b, Table 4c and 
Figure 3. 
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Table 3(a). Summary of Thermus oshimai JL-2 genome: one chromosome and two plasmids 



Label 


Size (Mb) 


Topology 


INSDC identifier 


RefSeq ID 


Chromosome 


2.072393 


Circular 


CP003249.1 




Plasmid pTHEOSOl 


0.271713 


Circular 


CP003250.1 




Plasmid pTHEOS02 


0.057223 


Circular 


CP003251.1 





Table 3(b). Nucleotide content and gene count levels of Thermus oshimai JL-2 genome 



Attribute 


Value 


% of Total 3 


Genome size (bp) 


2,401,329 


100.00 


DNA coding region (bp) 


2,251,025 


93.74 


DNA G+C content (bp) 


1,646,250 


68.56 


Total genes 


2,548 


100.00 


RNA genes 


60 


2.35 


Protein-coding genes 


2,488 


97.65 


Pseudogenes 


53 


2.08 


Genes in paralog clusters 


1,099 


43.13 


Genes with function prediction 


2,014 


79.04 


Genes assigned to COGs 


2,003 


78.61 


Genes assigned Pfam domains 


1,998 


78.41 


Genes with signal peptides 


862 


33.83 


Genes with transmembrane helices 


511 


20.05 


CRISPR repeats 


5 





a The total is based on either the size of the genome in base pairs or the total 
number of protein coding genes in the annotated genome. 

b Pseudogenes may also be counted as protein coding or RNA genes, so is 
not additive under total gene count. 



Table 3(c). Number of Thermus oshimai JL-2 genes associated with the 25 general COG functional categories 
Code Value %age Description 



J 


146 


6.67 


Translation 


A 


4 


0.18 


RNA processing and modification 


K 


114 


5.21 


Transcription 


L 


117 


5.35 


Replication, recombination and repair 


B 


2 


0.09 


Chromatin structure and dynamics 


D 


35 


1.60 


Cell cycle control, mitosis and meiosis 


Y 


0 


0 


Nuclear structure 


V 


25 


1.14 


Defense mechanisms 


T 


76 


3.47 


Signal transduction mechanisms 


M 


90 


4.11 


Cell wall/membrane biogenesis 


N 


23 


1.05 


Cell motility 


Z 


1 


0.05 


Cytoskeleton 


W 


0 


0 


Extracellular structures 


u 


44 


2.01 


Intracellular trafficking and secretion 


o 


85 


3.88 


Posttranslational modification, protein turnover, chaperones 


c 


154 


7.04 


Energy production and conversion 


G 


132 


6.03 


Carbohydrate transport and metabolism 


E 


219 


10.01 


Amino acid transport and metabolism 


F 


74 


3.38 


Nucleotide transport and metabolism 


H 


126 


5.76 


Coenzyme transport and metabolism 


I 


89 


4.07 


Lipid transport and metabolism 


P 


99 


4.52 


Inorganic ion transport and metabolism 


Q 


51 


2.33 


Secondary metabolites biosynthesis, transport and catabolism 


R 


289 


13.21 


General function prediction only 


S 


193 


8.82 


Function unknown 




545 


21.39 


Not in COGs 



'The total is based on the total number of protein coding genes in the annotated genome. 
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Thermus oshimai JL-2 

chromosome 500 kbp "I 
L isookbp (2, 072, 393 bp) 




RNA processing and modification (A) 

Chromatin structure and dynamics (B) 

Translation (J) 

Transcription (K) 

Replication and repair (L) 

Cell cycle control and mitosis (D) 

Post-translation modification, etc (O) 

Cell wall/membrana/envelop biogenesis (M) 

Cell motility (N) 

Inorganic Ion transport and metabolism (P) 

Signal Transduction (T) 

Intracellular trafficking, etc (U) 

Defense mechanisms (V) 

Extracellular structures (VU) 

Nuclear structure (Y) 

Cytoskeleton (Z) 

Energy production and conservation (C) 
Carbohydrate transport and metabolism (G) 
Coenzyme transport and metabolism (£) 



Lipid transport and metabolism (I) 
a Secondary metabolites, etc (O) 

General function prediction only (R) 

Function unknown IS) 
a No COGS identified 

CDS 

tRNA 

rRNA 

Other 

■ „ 100 % 

■ £ 98 % 

a 96% 

« 94 % 
3 92 % 
■ 90 % 
B 88 % 

86 % 
84 % 
82 % 
0 % 



Thermus thermophilus HB8 
: Thermus thermophilus HB27 
: Thermus scotoductus SA-01 



Figure 2. Map of T. oshimai JL-2 chromosome compared with other Thermus chromosomes. The outer four circles show the 
genes in forward and reverse strands and their corresponding COG categories. BLASTN hits (percentage identities) from T. 
thermophilus HB8 (1), T. thermophilus HB27 (2), and T. scotoductus SA-01 (3) chromosomes are shown in the inner three cir- 
cles. Maps were created using CGView Comparison Tool [32]. 



Table 4a. Summary of Thermus thermophilus JL-18 genome: one chromosome and two plasmids 



Label 


Size (Mb) 


Topology 


INSDC identifier 


RefSeq ID 


Chromosome 


1.902595 


Circular 


CP003252.1 


NC_01 7587.1 


Plasmid pTTJL1801 


0.265886 


Circular 


CP003253.1 


NC_01 7588.1 


Plasmid pTTJL1802 


0.0142731 


Circular 


CP003254.1 


NC_01 7590.1 



Comparison with other sequenced genomes 

The chromosome of T. thermophilus JL-18 was 
compared with the chromosomes of T. 
thermophilus strains HB8 and HB27 [38] using 
nucmer [39]. The megaplasmid pTTJL1801 was 
also compared with the megaplasmid sequences 
of HB8 and HB27. Dot plot results from this analy- 
sis [Figure 4(a]] demonstrate a high degree of 
synteny between the chromosomes of JL-18, HB8, 
and HB27; however, little synteny exists between 
the megaplasmids. T. oshimai JL-2 chromosome 
and megaplasmid sequences were also compared 



with those of T. thermophilus JL-18; however, little 
very synteny was apparent (Figure 4(b]]. 

Profiles of metabolic networks and pathways 

T. oshimai JL-2 and T. thermophilus JL-18 genomes 
encode genes for complete glycolysis, tricarboxylic 
acid (TCA] cycle, and pentose phosphate pathway 
(Figure 5). The genomes also encode glucosidases, 
glycosidases, proteases, and peptidases, highlight- 
ing the ability of these species to use various car- 
bohydrate and peptide substrates. Thus, central 
carbon metabolic pathways are very similar to 
those of T. thermophilus HB27 [38] and T. 
scotoductus SA-01 [41]. 
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Table 4b. Nucleotide content and gene count levels of Thermus thermophilus JL-18 genome 



Attribute 


Value 


% of total 3 


Genome size (bp) 


2,31 1,212 


100.00 


DNA coding region (bp) 


2,172,588 


94.00 


DNA G+C content (bp) 


1,594,227 


68.98 


Total genes b 


2,508 


100.00 


RNA genes 


56 


2.23 


Protein-coding genes 


2,452 


97.77 


Pseudogenes 


50 


1 .99 


Genes in paralog clusters 


1,069 


42.62 


Genes with function prediction 


1,979 


78.91 


Genes assigned to COGs 


1,992 


79.43 


Genes assigned Pfam domains 


1,962 


78.23 


Genes with signal peptides 


464 


18.5 


Genes with transmembrane helices 


518 


20.65 


CRISPR repeats 


3 




a The total is based on either the size of the genome in 
genes in the annotated genome. 

b Pseudogenes may also be counted as protein codin 


base pairs or the total 
g or RNA genes, so is 


number of protein coding 
not additive under total 



gene count. 




Thermus thermophilus JL- 
chromosome 
(1, 902, 595 bp) 




RNA processing and modification (A) 
Chromatin structure and dynamics (B) 
Translation (J) 



Call • 



.in,! f 



t (D) 



ion modification. etc (O) 
(I'll vi . ill mrmlir.mc envelop biogenesis (Ml 

Cell motility (N) 

Inorganic ion transport and metabolism (P) 

Signal Transduction ■: i i 

Intracellular trafficking, etc (U) 

Drlrnsr mechanisms (V) 

Extracellular structures (W) 

Nuclear structure (V) 

Cytoskeleton (Z) 

Energy production and conservation (CI 
Carbohydrate transport and metabolism <C) 
Coenzyme transport and metabolism (El 



Lipid transport and metabolism (I) 
' Secondary metabolites, etc (01 
General function prediction only (R) 



No COGs identified 



1: Thermus thermophilus HB8 
2: Thermus thermophilus HB27 



3: Thermus scotoductus SA-01 



Figure 3. Map of T. thermophilus JL-18 chromosome compared with other Thermus chromosomes. The outer 
four circles show the genes in forward and reverse strands and their corresponding COG categories. BLASTN hits 
(percentage identities) from T. thermophilus HB8 (1), T. thermophilus HB27 (2), and T. scotoductus SA-01 (3) 
chromosomes are shown in the inner three circles. Maps were created using CGView Comparison Tool [32]. 
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Code 


Value 


%age a 


Description 


J 


148 


6.79 


Translation 


A 


1 


0.05 


RNA processing and modification 


K 


104 


4.77 


Transcription 


L 


130 


5.97 


Replication, recombination and repair 


B 


2 


0.09 


Chromatin structure and dynamics 


D 


33 


1.51 


Cell cycle control, mitosis and meiosis 


Y 


0 


0 


Nuclear structure 


V 


25 


1.15 


Defense mechanisms 


T 


67 


3.07 


Signal transduction mechanisms 


M 


87 


3.99 


Cell wall/membrane biogenesis 


N 


30 


1.38 


Cell motility 


Z 


1 


0.05 


Cytoskeleton 


W 


0 


0 


Extracellular structures 


u 


57 


2.62 


Intracellular trafficking and secretion 


o 


82 


3.76 


Posttranslational modification, protein turnover, chaperones 


c 


149 


6.84 


Energy production and conversion 


G 


125 


5.74 


Carbohydrate transport and metabolism 


E 


216 


9.91 


Amino acid transport and metabolism 


F 


64 


2.94 


Nucleotide transport and metabolism 


H 


119 


5.46 


Coenzyme transport and metabolism 


I 


94 


4.31 


Lipid transport and metabolism 


P 


96 


4.41 


Inorganic ion transport and metabolism 


Q 


57 


2.62 


Secondary metabolites biosynthesis, transport and catabolism 


R 


291 


13.35 


General function prediction only 


S 


201 


9.22 


Function unknown 




516 


20.57 


Not in COGs 



a The total is based on the total number of protein coding genes in the annotated genome. 



Genes involved in denitrification 

Denitrification involves the conversion of nitrate 
to dinitrogen through the intermediates nitrite, 
nitric oxide, and nitrous oxide and is mediated by 
nar, nir, nor, and nos genes [4]. Incomplete 
denitrification phenotypes terminating in the pro- 
duction of nitrous oxide have recently been re- 
ported for a large number of Thermus isolates, in- 
cluding T. oshimai JL-2 and T. thermophilus JL-18 
[6]. 

Figure 6 shows the organization of the nar operon 
and neighboring genes involved in denitrification 
in T. oshimai JL-2, T. thermophilus JL-18, and T. 
scotoductus SA-01. These gene clusters are located 
on the megaplasmids of T. oshimai JL-2 and T. 



thermophilus JL-18, as in other T. thermophilus 
strains [44,45]. They are located on the chromo- 
some in T. scotoductus SA-01 [41]. The nar oper- 
ons show a high degree of synteny and all include 
genes encoding the membrane-bound nitrate 
reductase (NarGHI], the associated periplasmic 
cytochrome NarC, and the dedicated chaperone 
NarJ. All three strains contained homologs of 
NarKl, which is a member of the major facilitator 
superfamily that likely functions as a ni- 
trate/proton symporter [46,47]. However, some 
experiments in T. thermophilus HB8 suggest 
NarKl might also function in nitrite extrusion 
[39]. T. oshimai JL-2 and T. scotoductus SA-01 also 
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contain homologs of NarK2 (annotated as nep in T. 
scotoductus SA-01 [41]], which likely encodes a 
nitrate/nitrite antiporter [44,48]. No significant 
BLASTP hits for periplasmic nitrate reductase 
subunits NapB and NapC were found in T. oshimai 
JL-2 and T. thermophilus JL-18, consistent with the 
use of the Nar system in the Thermales. 

All three strains contain a dnrST operon adjacent 
to, but divergently transcribed from, the narGHJIK 
operon. dnrST encodes transcriptional activators 
responsible for upregulation of the nitrate respira- 
tion pathway in the absence of O2 and the pres- 
ence of nitrogen oxides or oxyanions [42] (Figure 
6). 

Both the species contain a putative nirK, which 
encodes the NO-forming, Cu-containing nitrite 
reductase. In addition, T. oshimai JL-2 and T. 
scotoductus SA-01 both harbor nirS [41], which 
encodes the isofunctional tetraheme cytochrome 



cdi-containing nitrite reductase. Previous studies 
have suggested that bacteria use either NirK or 
NirS, but not both, for the reduction of nitrite [49]. 
The unique presence of NirK and NirS in T. 
oshimai JL-2 and T. scotoductus SA-01 likely en- 
hances their denitrification abilities since 
isoenzymes are typically kinetically distinct 
and/or regulated differently. This idea is con- 
sistent with the distinct denitrification pheno- 
types of T. oshimai strains as compared to T. 
thermophilus strains reported previously, includ- 
ing strains T. oshimai JL-2 and T. thermophilus JL- 
18 [6]. In those studies, nitrite accumulated in the 
medium at concentrations of <150 |iM in T. 
thermophilus strains, whereas it was rapidly pro- 
duced to concentrations >200 |iM but consumed 
rapidly to below method detection limits in T. 
oshimai strains. 




Thermits thermophilus HB8 Thermus thermophilus HB27 




Thermus thermophilus HB8 Thermus thermophilus HB27 

plasmid pTT27 plasmid pTT27 



Figure 4(a). Dot plot comparison of T. thermophilus JL-18 chromosome and megaplasmid DNA sequence with 
those of the strains HB8 and HB27. 
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Figure 4(b). Dot plot comparing the chromosome and megaplasmid DNA sequence of T. oshimai JL-2 and T. 
thermophilus JL-18. 



NirK functions as a homo-trimer [50] and contains 
type 1 (blue] and type 2 (non-blue] copper-binding 
residues [49]. Comparison of the NirK from T. 
oshimai JL-2 and T. scotoductus SA-01 with previ- 
ously studied NirK amino acid sequences revealed 
that six of the seven copper-binding residues are 
conserved, except for a single methionine (M] to 
glutamine (Q] substitution in both Thermus pro- 
teins (Figure 7; indicated by an asterisk (*]]. Glu- 
tamine, not methionine, is the copper-binding lig- 
and in the case of stellacyanin, a blue (type 1] cop- 
per-containing protein [52,53]. A M121Q recombi- 
nant protein of Alcaligenes denitrificans azurin 
showed similar electron paramagnetic resonance 
(EPR], but exhibited a 100-fold lower redox activity 
when compared to wild-type azurin [54]. There- 
fore, although the methionine is replaced with a 
glutamine in the T. oshimai JL-2 NirK, it is possible 
that this glutamine residue can function as a cop- 
per-binding ligand similar to stellacyanin and 
azurin. The large and small subunits of nitric oxide 
reductase (NorB and NorC] are predicted to be 



co-transcribed along with nitrite reductases in T 
oshimai JL-2, T. thermophilus JL-18 and T. 
scotoductus SA-01 (Figure 6]. 

Genes encoding the 15 subunit NADH-quinone 
oxidoreductase [55] were identified in both ge- 
nomes (Theos_0703 to 0716, 1811 in T. oshimai JL- 
2; TTJL18JL786 to 1799, 1580 T. thermophilus JL- 
18]. nrcDEFN, a four gene operon encoding a novel 
NADH dehydrogenase, is adjacent to the nar oper- 
on in the megaplasmid of T. thermophilus HB8 and 
has been previously implicated in nitrate reduction 
[43]. In T. thermophilus JL-18, the operon is present 
(Figure 6], although (TTJL18_2313] is truncated 
(NarE in HB8: 232 AA, in JL-18: 78 AAJ. In T. 
oshimai JL-2, only nrcN is present. Theos_0161 and 
Theos_0162, orthologs of Wolinella succinogenes 
NrfA and NrfH [56], respectively, were identified in 
T. oshimai JL-2 suggesting that T. oshimai JL-2 may 
be capable of respiratory nitrite ammonification, 
although this phenotype has not yet been observed 
in Thermus [6]. 
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Other possible electron transport components in- 
clude a ftcz3-type heme-copper oxidase 
(Theos_1499, 1498, 1497, T. oshimai JL-2; 
TTJL18_0925, 0926, 0927 T. thermophilus JL-18) 
and bc\ complex encoded by the FbcCDFB operon 
[57]. (Theos_0106 to 0109, T. oshimai JL-2; 
TTJL18_2018 to 2021 T. thermophilus JL-18). In 
addition, both T. oshimai JL-2 and T. thermophilus 
JL-18 harbor genes for archaeal-type Vo-Vi (vacuo- 
lar) type ATPases, which appears to have been 
acquired from Archaea prior to the divergence of 
the modern Thermales [58]. 

Genes involved in iron reduction 

T. scotoductus SA-01 has been reported to be ca- 
pable of dissimilatory Fe 3+ reduction; however, 
the biochemical basis of iron reduction has not 
been elucidated in Thermus [41,59]. Sequences of 
proteins involved in iron reduction [60] in 
Shewanella oneidensis MR-1 (MtrA, MtrF, OmcA) 
and Geobacter sulfurreducens KN400 (OmcB, 
OmcE, OmcS, OmcT, OmcZ) were used as search 



queries into Thermus genomes using BLASTP. No 
hits were found in T. oshimai JL-2, T. thermophilus 
JL-18, or T. scotoductus SA-01. This suggests that 
the biochemical basis of iron reduction is distinct 
in Thermus compared to Shewanella and 
Geobacter, and offers no predictive information on 
whether T. oshimai JL-2 and T. thermophilus JL-18 
may be able to respire iron. 

Genes involved in sulfur oxidation 

A complete sox cluster comprising of 15 genes, 
including soxCD, is present in T oshimai JL-2 and 
T. thermophilus JL-18 genomes. SoxCD is essential 
for chemotrophic growth of P. pantotrophus [61]. 
Taken together, this suggests that T. oshimai JL-2 
and T. thermophilus JL-18 may use thiosulfate as 
an electron donor and are similar to other sulfur- 
oxidizing Thermus strains including T. scotoductus 
IT-7254 [62] and T. scotoductus SA-01 [41]. Other 
T. thermophilus genomes also harbor this gene 
cluster, suggesting thiosulfate oxidation may be 
widely distributed in Thermus [38]. 




Figure 5. Metabolic pathways identified using iPATH2 [40]. Orange lines are common pathways that were identified in T. 
oshimai JL-2 and J. thermophilus JL-18. Blue lines indicate pathways unique to T. oshimai JL-2 and red lines indicate pathways 
unique to T. thermophilus JL-18. 
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Figure 6. Map showing the organization of nar operon and neighboring genes involved in denitrification located on the 
megaplasmids of T. oshimai JL-2 (pTHEOSOl) and T. thermophilus\LAQ (pTT|L180l) and the chromosome of T. 
scotoductus SA-01. Fe: heme protein-containing nitrite reductase, Cu: copper-containing nitrite reductase. Numbers be- 
low the genes indicate the provisional ORF numbers in T. oshimai JL-2 (Theos_1057 - Theos_1036) and T. thermophilus 
JL-18 (TtJL18_2297 to TtJL1 8_2327), the locations in the megaplasmid are indicated below, nar. nitrate reductase; nir. 
nitrite reductase; nos: nitric oxidereductase; dnr. denitrification regulator [41-43]. 
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P. aureofaciens 
A. cycloclastes 
R. sphaeroidcs 1 7025 
R. sphaeroides KD131 
A. faecalis 
A. xylosoxidans 
N. europaea 
N. meningitidis Z2491 
T. oshimai JL-2 
T. scotoductus SA-01 

Figure 7. Thermus oshimai JL-2 gene Theos_1053 encodes a Copper-containing nitrite reductase. Amino acid sequences 
of known Cu-containing nitrite reductases from Pseudomonas aureofaciens (P. aureofaciens, Gl: 287907), Achromobacter 
cycloclastes (A. cycloclastes, Gl: 157835402), Rhodobacter sphaeroides ATCC 17025 (R. sphaeroides 17025, Gl: 
146277634), Rhodobacter sphaeroides KD131 (R. sphaeroides KD131, Gl: 221638756), Alcaligenes faecalis (A. faecalis, 
Gl: 393758960), Alcaligenes xylosoxidans {A. xylosoxidans, Gl: 422318032), Nitrosomonas europaea (N. europaea, Gl: 
30248928), Neisseria meningitidis Z2491 (N. meningitidis Z2491, Gl: 218768658) and Thermus scotoductus SA-01 (T. 
scotoductus SA-01, Gl: 320450829) were aligned using Muscle v3.8.31 [51] along with Thermus oshimai JL-2 (7. oshimai 
JL-2, Gl: 410732282) Theos_1053. Putative copper-binding residues are indicated with downward arrows according to 
their classes: 1 : type 1 (blue) Cu; 2: type 2 (nonblue) Cu [49]. Numbers on left and right of the alignments refer to positions 
in the alignment. Asterisk (*) indicates the MfJQ substitution ifi oshimai JL-2 and T. scotoductus SA-01 . 
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A variety of chemotrophs and anoxygenic 
phototrophs can oxidize hydrogen sulfide, organic 
sulfur compounds, sulfite, and thiosulfate as elec- 
tron donors for respiration [63]. Reconstituted pro- 
teins of SoxXA, SoxYZ, SoxB and SoxCD together, 
but not alone, mediate the oxidation of thiosulfate, 
sulfite, sulfur, and hydrogen sulfide in Paratrophus 
pantotrophus [61]. The absence of free intermedi- 
ates of sulfur oxidation and the occurrence of sul- 
fite oxidation without SoxCD in P. pantotrophus ex- 
cludes SoxCD as a sulfite dehydrogenase and pro- 
vides evidence to its role as a sulfur dehydrogenase 
with protein-bound sulfur atom [61]. 

Polysulfide reductase in T. oshimai JL-2 

In T. oshimai JL-2, three proteins showed high se- 
quence identity to PsrA (88%; Theos_0751), PsrB 
(86%; Theos_0750], and PsrC (83%; Theos_0749) 
of T. thermophilus HB27, which is likely involved 
in anaerobic respiration using polysulfide as a 
terminal electron acceptor. In T. thermophilus 
HB27, PsrA is the putative catalytic subunit con- 
taining two molybdopterin guanine dinucleotide 
co-factors and a cubane-type [4Fe-4S] cluster. 
Electron transfer is likely mediated by PsrB, which 
also contains a [4Fe-4S] cluster, while PsrC is a 
putative transmembrane protein that contains the 
electron carrier menaquinone-7 (MK-7]. PSR func- 
tions as a hexamer (composed of 2 subunits each 
of A, B and C] and catalyzes the reactions: 
MKH 2 ->MK + 2H + + 2e- in the membrane, and S n 2 -+ 
2e- + 2H + + Sn i 2 " + H2S in the periplasm [64]. How- 
ever, the Thermus PsrABC proteins exhibit very 
low identity to Wolinella succinogenes PsrABC pro- 
teins that have been functionally characterized 
(PsrA: 33%, PsrB 46%, no clear BLASTP hits found 
in T. oshimai JL-2 for W. succinogenes PsrC) [65]. 
In Wolinella succinogenes, formate dehydrogenase 
or hydrogenase and polysulfide reductase form 
the electron transport chain and mediate the re- 
duction of polysulfide with formate or H2 [64]. In 
T. oshimai JL-2, Theos_1377 encodes a putative 
formate dehydrogenase alpha subunit. Another 
gene, Theos_llll, encodes a putative formate de- 
hydrogenase family accessory protein (FdhD], 
which is required for regulation of the formate 
dehydrogenase catalytic subunit [66] and is con- 
served in many members of the Thermaceae, in- 
cluding T. scotoductus SA-01 (TSC_cl0040]. Alt- 
hough the genes needed for polysulfide reduction 
are present, polysulfide reduction in T. oshimai JL- 
2 has not been tested. 



Genes involved in DNA uptake 

A significant number of genes in hyperthermophilic 
bacteria are of archaeal origin, and appear to have 
been acquired through inter-domain gene transfer 
[67], which is mediated by both transformation and 
conjugation systems [68]. T. thermophilus HB27 is 
naturally competent to both linear and circular 
DNA, and DNA transport mechanisms in this spe- 
cies have been well studied [69,70]. The genome of 
T. oshimai JL-2 and T. thermophilus JL-18 both con- 
tain homologs of DNA transport genes (Table 5], 
suggesting that both T. oshimai JL-2 and T. 
thermophilus JL-18 are naturally competent. 

Conclusions 

We report the finished genomes of T. oshimai JL-2 
and T. thermophilus JL-18. T. oshimai JL-2 is the first 
complete genome to be reported for this species, 
while T. thermophilus JL-18 is the fourth genome to 
be reported for T. thermophilus. Analysis of the ge- 
nomes revealed that they encode enzymes for the 
reduction of nitrate to nitrous oxide, which is con- 
sistent with the high flux of nitrous oxide reported 
in GBS [6], and explains the truncated 
denitrification phenotype reported for many 
Thermus isolates obtained from that system [6]. It 
is intriguing that Thermus scotoductus SA-01 also 
has genes encoding the sequential reduction of ni- 
trate to nitrous oxide but lacks genes encoding the 
nitrous oxide reductase. The high degree of synteny 
in the respiratory gene cluster combined with the 
conserved absence of the nitrous oxide reductase 
suggests incomplete denitrification might be a pre- 
viously unrecognized but conserved feature of 
denitrification pathways in the genus Thermus, alt- 
hough T. thermophilus NAR1 appears to be capable 
of complete denitrification to N2 [73]. Another unu- 
sual feature of the T. oshimai JL-2 and T. 
scotoductus SA-01 denitrification systems is the 
apparent presence of the NO-forming, Cu- 
containing nitrite reductase, NirK, and the 
isofunctional tetraheme cytochrome cdi-containing 
nitrite reductase, NirS. 

T. oshimai JL-2 and T. thermophilus JL-18 also may 
be capable of sulfur oxidation since they both en- 
code a complete, chromosomal sox cluster. Howev- 
er, experiments with GBS sediments failed to 
demonstrate a stimulation of denitrification when 
thiosulfate was added in excess [74], suggesting 
thiosulfate oxidation may not be coupled to 
denitrification in these organisms. The presence of 
psrA, psrB and psrC genes encoding polysulfide 
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reducatase in T. oshimai JL-2 suggests the ability to 
reduce polysulfide. The function of these putative 
pathways could be tested with pure cultures in the 
laboratory. 

The presence of complete macromolecular ma- 
chinery for natural competence and the presence 
of megaplasmids harboring genes for ni- 
trate/nitrite reduction and thermophily points out 
that T. oshimai JL-2 and T. thermophilus JL-18 
could have acquired innumerable genes through 



intra- and inter-domain gene transfer, and sug- 
gests considerable plasticity in denitrification 
pathways. Considering the importance of these 
organisms in the nitrogen biogeochemical cycle, 
and their potential as sources of enzymes for bio- 
technology applications, the complete genome se- 
quences of T. oshimai JL-2 and T. thermophilus JL- 
18 are valuable resources for both basic and ap- 
plied research- 



Table 5. Identification of competence proteins in T. oshimai JL-2 and T. thermophilus JL-1 8 by IMG/ER [71 ] . + 
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Linkage of periplasmic and cytoplasmic proteins 


PilQ 


Theos_0435 


TtJLI f 


3_0665 


Directing DNA transporter through OM 


ComZ 


Theos_1239 


TtJLI { 


3_0832 


IM protein, function unknown 


PilM 


Theos_0439 


TtJLI { 


3_0669 


ATPase, function unknown 


PilN 


Theos_0438 


TtJLI { 


3_0668 


IM protein, function unknown 


PilO 


Theos_0437 


TtJLI { 


3_0667 


IM protein, function unknown 


PilW 


Theos_0436 


TtJLI { 


3_0666 


OM protein, stabilization of PilQ 



+ BLASTP analysis using sequences 
modified from [72]. 



of known competence proteins from T. thermophilus HB27 as queries. Table 
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